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ABSTRACT 

We map the stellar structure of the Galactic thick disk and halo by applying color-magnitude 
diagram (CMD) fitting to photometric data from the SEGUE survey. The SEGUE imaging scans 
allow, for the first time, a comprehensive analysis of Milky Way structure at both high and low latitudes 
using uniform SDSS photometry. Incorporating photometry of all relevant stars simultaneously, CMD 
fitting bypasses the need to choose single tracer populations. Using old stellar populations of differing 
metallicities as templates we obtain a sparse three-dimensional map of the stellar mass distribution at 
\Z\ >1 kpc. Fitting a smooth Milky Way model comprising exponential thin and thick disks and an 
axisymmetric power-law halo allows us to constrain the structural parameters of the thick disk and 
halo. The thick-disk scale height and length are well constrained at 0.75 ± 0.07 kpc and 4.1 ± 0.4 kpc, 
respectively. We find a stellar halo flattening within ^25 kpc of c/a = 0.88 ± 0.03 and a power-law 
index of 2.75 ± 0.07 (for 7 kpc < Rgc ^ 30 kpc). The model fits yield thick-disk and stellar halo 
densities at the solar location of p t Mck,& = 10~ 2 - 3±01 M Q pc" 3 and phaio,& = 10~ 4 - 20±0 05 M pc -3 , 
averaging over any substructures. Our analysis provides the first clear in situ evidence for a radial 
metallicity gradient in the Milky Way's stellar halo: within R <15 kpc the stellar halo has a mean 
metallicity of [Fe/H] ~ —1 .6, which shifts to [F e/H] ~ —2.2 at larger radii, in line with the two- 
component halo deduced bv lCarollo et al.l (|2007|) from a local kinematic analysis. Subtraction of the 
best-fit smooth and symmetric model from the overall density maps reveals a wealth of substructures 
at all latitudes, some attributable to known streams and overdensities, and some new. A simple warp 
cannot account for the low latitude substructure, as overdensities occur simultaneously above and 
below the Galactic plane. 

Subject headings: Galaxy: structure - Galaxy: disk - Galaxy: halo - stars: statistics 



1. INTRODUCTION 

Studying the stellar structure of the Milky Way has a 
long history, as early models based on star counts were 
already constru cted by the likes of iHerschel (|1785[ ) and 
iKaptevnl (|1922f ). Today we know that four 'components' 
make for a sensible approximate description of the Milky 
Way stellar body, each with different structural, kine- 
matic and population characteristics: the bulge, the thin 
and the thick disk, and the stellar halo. Accurate de- 
termination of the properties of these components is a 
difficult undertaking, as it requires surveys with (1) suf- 
ficient sky coverage to assess the overall geometry, (2) 
sufficient depth for mapping stars to ~10 to 30 kpc, and 
(3) sufficient information (e.g., color-based luminosity es- 
timates) to obtain reasonable distance estimates for these 
stars. Furthermore, the presence of dust in the plane of 
the disk clouds our view of most of the Galaxy at wave- 
lengths shorter than a few microns. For this reason, most 
previous efforts to study Milky Way structure at low lat- 
itu des have us e d (ne ar-) infrared data, such as the work 
by Kent et all (|1991l) based on data from the Spacelab 
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infrared t elescope, or more recent work using 2MASS 
data (e.g. iMomanv etaTl 120061: iRevle et al.ll2009f > or the 
GLIMPSE survey (jBeniamin et al.ll2005| ) While near- 
infrared surveys provide the only way of studying the low 
latitude regions, because they suffer much less from dust 
extinction, they are currently less sensitive than their 
optical counterparts and often rely on tracer populations 
with less accurate distances. 

Over the past f ew years, data from 2MASS 
(ISkrutskie et al.1 12006D and the Sloan Digital Sky Sur- 
vev (SDSS: lYork et al.l[2000HAbazaiian et alJl2009h have 
played a key role in revolutionizing our empirical picture 
of the stellar components of the Milky Way. As the SDSS 
survey avoided the Galactic plane, much of the work us- 
ing SDSS data has focused away from the bulk of the 
Milk y Way's stars to h igh galactic latitudes. Nonethe 
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less, I Juric et "aLl (|2008f ) have been able to constrain the 
structural parameters of the stellar halo and the thin 
and thick disk. Beyond the global parameters, these sur- 
veys have also led to the discovery of substructure in 
both the stellar disk dNewberg et al.l 120021: llbata et al.l 
120031 iMartin et al.ll2q04t iBellazzini et al ll2006l) and halo 
dNewberg et al.l l2002t jBelqkurov et al.1 I^OOajGrillmairl 
| 2006t iBelokurov et al.H2007t Uuric et alJl2008t feell et al.1 
2008). This detailed substructure provides information 
on how the Milky Way formed and evolved; as galax- 
ies assemble through mergers and accretions of smaller 
systems, these events leave identifiable signatures in their 
structure and kinematics. This also has cosmological im- 
plications, as the structure of the dark halo and the in- 
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teractions with dark subhalos c an significantly influence 
the appearance of a galaxy (e.g. iKazantzidis et alJl2008t 
lYounger et alll2008l) . 

In this paper we derive a new map of the stellar dis- 
tribution in the Milky Way, in order to study the large- 
scale structure of the thick disk and stellar halo. Though 
covering only the Northern celestial hemisphere, and 
that only sparsely, this map covers both high and low 
latitudes and uses rigorously derived stellar distances. 
This map is based on one of the constituent projects 
of the extended SDSS survey (SDSS II), SEGUE (Sloan 
Extension for Gala ctic Understanding and Exploration; 
lYannv et al1l2009aft . which is an imaging and spectro- 
scopic survey aimed at the study of the Milky Way and 
its stellar populations. Whereas the main SDSS survey 
avoided low Galactic latitudes, SEGUE imaging scans 
go through the Galactic plane (see Figure [TJ), allowing 
a 'picket- fence' view of the Galaxy at these low lati- 
tudes. It is this deep, uniform data set that enables 
a systematic analysis of stellar structure as function of 
Galactic latitude and longitude, which we base on color- 
m agnitude diagram (CM P) fitting techniques developed 
bv lde Jong et al.1 ([2008). A crucial advantage of an anal- 
ysis based on CMDs is that they bypass the often ad- 
hoc and sometimes problematic choices of stellar 'tracers' 
(such as K-dwarfs or F-stars) to delineate the structure. 
Instead, this procedure results immediately in a mass- 
density map for stellar population components of, say, 
different metallicity. 

To avoid having to fit the complex stellar population 
make-up of the thin disk, we in the end consider only 
stars that are more than 1 kpc away from the Galactic 
mid-plane, and treat the thin disk as a small perturba- 
tion to our results. As SEGUE is a northern hemisphere 
survey, it stays far away from the Galactic bulge, which 
can be safely ignored. On the other hand, the resulting 
maps of stellar mass density in the thick disk and halo 
allow us to measure the structural parameters of these 
components. In addition, we show that subtracting a 
smooth Galactic model reveals a wealth of substructure, 
especially at low latitudes, which we will study in detail 
in a forthcoming publication. 

The outline of this paper is as follows. Section [2] de- 
scribes the data used for our analysis, and Sj3] describes 
the CMD-fitting techniques we utilize. The resulting stel- 
lar distributions are presented in 21 where we also fit 
smooth Galactic models to them, derive structural pa- 
rameters, and look at deviations from the best fitting 
model. Finally, we summarize and discuss our results in 



2. SEGUE PHOTOMETRY AND COVERAGE 

As part of the SEGUE survey (|Yannv et al.ll2009aD . 
ten 2.5°-wide scans, crossing the Galactic plane at fixed 
longit udes, were imaged in five passbands (u, g, r, i, 
and 2: lGurm^t^fl99al200llHogg et al.ll2001h . These 
data are available as part of SDSS Data Release 7 
(DR7, lAbazaiian et al.l [2009ft . An automated data re- 
duction and analysis pipeline produces accurate astro- 
metric and photometric measurements for all detected 
objec t s dLupton. Gunn. fc Szalav 1999; Stoughto n et al.l 
2002t iSmith et all 120021: iPier et alTl2003t llvezic et al.1 
2004 iTucker et al.l 120061 ) ■ reaching a photo metric accu- 
racy of 2% down to g ~22.5 and r -22 (I vezic et al.l 
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Fig. 1. — Overvie w of low-latitude imaging scans from SEGUE 
llYannv et al.ll2009al ) used for this analysis, shown as the vertical 
lines in Galactic coordinates. The grayscale shows dust red dening 
estimates based on the extinction maps from Schlcgcl ct al. (1998), 
with increasingly dark regions corresponding to E(B — V) >0.1, 
0.25, 0.5 and 1.0 mag. Only regions with E(B — V) < 0.25 are 
used in the analysis. 



[2001 iSesar et all [2007T ) . Figure Q] shows the coverage of 
the scans between Galactic latitudes of +50° and —50°. 
Most scans have some overlap with the SDSS main sur- 
vey or SDSS Legacy area in the North Galactic cap, and 
we extract 2.5°-wide strips of data from DR7 to extend 
the photometry to Galactic latitudes of 85°. For the 
CMD fitting analysis we restrict ourselves to two bands, 
g and r. These two bands are the most useful for our 
population/distance analysis because, among the three 
photometrically most sensitive bands (<?, r and i), they 
combine to give the largest offset in magnitude between 
stellar main sequences of different metallicities. The 
SEGUE imaging in DR7 have been photometrically cal- 
ibrate d following the procedure of iPadmanabhan et al.l 
(2008). The procedure basically uses the overlap of the 
SEGUE data scans with the Legacy SDSS imaging scans 
to solve for the photometric zero-points and flat-field co- 
efficients of the SDSS 2.5m camera system. This overlap 
generally occurs in uncrowded regions of sky (|6| > 20°) 
even for scans which extend to cross the Galactic plane. 
The photometric solution for the SEGUE scans was done 
incrementally rather than globally after the global solu- 
tion for all the SDSS Legacy scans. The accuracy of 
the photometry is estimated to be to about 1-5% in 
and r (< 2% in color). See §2.2 of lYannv etafl ((2009. 
for more details on the SEGUE imaging and calibra- 
tion. We de-redden all data using the dust maps from 
ISchlegel et al.l (119981). incl uding the correction suggested 
bv iBonifacio et al.l <|2000T >. provided through the SDSS 
catalog server. As most of the dust is confined close to 
the Galactic plane, it is in the foreground for all stars 
more distant than 1 kpc from the mid-plane which we 
employ in our analysis. 

To illustrate the basis for our subsequent modeling, 
several de-reddened CMDs drawn from different latitudes 
in the SEGUE stripe at Z=94°are shown in Figured The 
most obvious difference between different latitudes is the 
variation in reddening, which ranges from an E(B — V) 
of up to several magnitudes very close to or in the plane 
of the Galaxy to a few hundredths of a magnitude at high 
latitudes. Once the reddening reaches several tenths of 
magnitudes (e.g., the upper right panel of Fig. [2]), the 
unaccounted differential reddening along the line of sight 
smears out the CMD features. For this reason, and be- 
cause the dereddening itself is less reliable in high extinc- 
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tion regions, we exclude regions with E(B — V) >0.25 
from our analysis; the excluded regions correspond to 
the three darkest shadings in Fig. [TJ Hence, the two 
lower panels in Figure [2] are representative for the CMDs 
used for this study. This cut in foreground extinction 
also means that we exclude the most crowded regions 
and thereby avoid possible problems due to crowding. 

Before carrying out a more formal CMD analysis, it 
is instructive to describe the main features seen in the 
CMDs and to realize which kind of stars they contain. 
Figure [3] shows the number density of stars in bins in the 
color-magnitude plane, or Hess diagram, of all stars in 
the SEGUE scan at Z=94° and b >30°. The large con- 
centration of stars in the lower right, i.e. at faint mag- 
nitudes and centered at a g — r color of ~1.4 consists of 
nearby, intrinsically faint and low-mass dwarf stars in the 
disk. The other plume of stars at g — r ~0.4 is where the 
main-sequence turn-off (MSTO) stars of old (>10 Gyr) 
stellar populations lie. A shift in the average color of this 
plume of stars is visible around g —18. This is caused 
by the difference in metallicity between the thick disk, 
which dominates at brighter magnitudes, and the Galac- 
tic halo, which dominates at larger distances. The color 
of the MSTO is a useful proxy for the metallicity of the 
dominant population at a given distance (e. g. see the dis- 
cussi on about this for SDSS isochrones in iGirardi et al.l 
2004). Assuming a regular and smooth Milky Way, the 
main change between different latitudes are the relative 
contributions of (thick) disk and halo stars to the CMD. 
For example, comparing the CMDs corresponding to lat- 
itudes of 30° and 70° (lower panels of Fig. [2]) shows 
how the number of disk stars decreases with respect to 
the halo stars. The field at 70° latitude adds complexity 
to this description, as it contains an example of stel- 
lar substructure; the main-sequence feature produced by 
the Sagittarius stream is readily visible, running from 
( ff -r, 5 )~(0.2,20.0) to (0.5,22.0). ' 

For our analysis of stellar structure we divide the 2.5°- 
wide SEGUE stripes into bins of 1° width in latitude from 
which CMDs are created. Each CMD is then analysed 
as described in the next Section. 



3. DISTANCE FITTING OF CMDS 

In most applications of CMD fitting, observed photom- 
etry is compared with models in order to determine which 
combination of simple stellar populations best resembles 
the data, thereby providing a model of the star forma- 
tion history (SFH) and age- metallicity relation (AMR). 
For the work presented here, we fit model stellar popula- 
tion CM Ds based on the iso chrones in SDSS filters pro- 
vided bv lGirardi et al.l (|2004l) to the SEGUE pho tometry 
The softwa re package we use for this (MATCH, iDolphinl 
[1993, HqoH) uses a maximum-likelihood method to deter- 
mine the best linear combination of models, after trans- 
forming both models and data to Hess diagrams (2D his- 
tograms of stars in the color- magnitude plane) , enabling 
a pixel-by-pixel comparison. For a proper comparison 
the models need to be convolved with a realistic photo- 
metric error and com pleteness model. The model we use 
here is the same as in lde Jong et al.1 (|2008f ) , based on the 
photo metric errors from the SDSS pipeline (jlvezic et al.l 
I2004T) and the completeness determined from a compari- 
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Fig. 2. — Typical color-magnitude diagrams taken from the 
SEGUE imaging s can at £=94°, deredde ned based on the dust ex- 
tinction maps from Schlcgcl ct al. ( 1998). In the upper right corner 
of each pa nel the Galactic lat itude and the average reddening in 
E(B — V) (Schlcgcl ct al. 1998) corresponding to each CMD are in- 
dicated. The latitude range from which the GMDs were extracted 
were chosen so that each contains approximately the same number 
of stars. 




Fig. 3. — A Hess diagram (greyscale) of the SEGUE scan 
at Z=94°and b >30°, in dereddened (g,g — r). Overplotted are 
isochrones corresponding to the three template stellar popula- 
tions fit to the data. From right to left the isochrones are for 
the 'thick-disk-like' ([Fe/H]=— 0.7), the 'inner-halo-like' intermedi- 
ate metallicity ([Fe/H]=— 1.3) and the 'outer-halo-like' metal-poor 
([Fe/H]=-2.2) populations. 

son 5 between SDSS and COMBO-17 6 . Since foreground 
reddening precludes the use of the most crowded regions 
in the mid-plane, crowding does not significantly affect 
the completeness in the data used here. 

Rather than determining the SFH and AMR, our aim 
in the present context is to map the structure of the 
stellar populations in the inner and outer Milky Way. 

5 See http://www.sdss.org/dr5/products/general/completeness.html 
e http:/ /www. mpia.de/COMBO/combo_index. html 
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(a) Observed Hess diagram (b) Model Hess diagram 
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(c) Residuals 
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Fig. 4. — Example of distance distribution fitting, based on the Hess/CMD diagram at 1=94° and f>=+30°; note that only the blue 
portion of the full color range in Figure[3]is shown here. No attempt to model low-mass thin disk stars dominating at g — r >0.8 is made 
here, (a) Observed distribution of stars in the color- magnitude plane, (b) Best-fit model to (a), consisting of a combination of template 
stellar populations at different distances, (c) Residuals after subtracting the model (b) from the data (a). Color bars indicate the number 
of stars in the Hess diagram bins. The mean residuals across the CMD are ~25% of the mean star counts across the entire CMD. 



For this, we use the dis tance-fitting option o f MATCH, 
described and tested in Ide Jong et al.1 ((2008) . To limit 
the number of free parameters and parameter degener- 
ation in the fits, we define a small set of template stel- 
lar populations. The choice of templates is motivated 
both by the populations we expect to be present in the 
outer disk and inner halo, and by the CMD features 
we have to fit. As discussed in iJU the most obvious 
indicator of population differences will be the MSTO 
color, which means the templates should span the ex- 
pected color range, preferably in regular intervals. The 
MSTO color of a stellar population depends both on its 
age and metallicity, but since both the thick disk and the 
halo stellar populations are known to be old, we choose 
to use a fixed age range and use metallicity offsets to 
probe the MSTO color range. We adopt three different 
metallicity bins, [Fe/H]=— 0.7, —1.3 and —2.2, respec- 
tively. All template populations have the same age range 
of 10.1<log[t/yrs]<10.2, and an assumed binary fraction 
of 0.5. For convenience, we choose a nomenclature for 
these components that is based on previous work: we re- 
fer to [F e/H]=— 0.7 as a thick -disk-like population and, 
following ICarollo et all (|2007ft . refer to the [Fe/H]=-1.3 
and —2.2 templates as the two halo populations, namely, 
an inner-halo- like and outer-halo-like population, respec- 
tively. Note, however, that in our subsequent analysis 
this choice of terms neither immediately prejudges the 
geometry, nor implies that these components are 'dis- 
tinct'. We simply presume that these three components 
are sufficient to completely describe the stellar popula- 
tions at any point that is at least 1 kpc above or be- 
low the Galactic mid-plane. Isochrones corresponding to 
these population templates are overplotted on Figure [3J 
The colors of the MSTO of populations 1 and 2 match 
the color of the prominent upper and lower plume of 
MSTO stars in Figure [3J respectively, while the third 
population has a bluer MSTO. Using a coarse grid of 
metallicities means that for stars with a metallicity dif- 
ferent from any of the template populations a discrepant 
distance is inferred, as the brightness of the MS depends 
is metallicity dependent (see Fig. [3]). The magnitude 
offsets in this case are ~0.3 and ~0.5 magnitudes, going 



from the highest ([Fe/H]=— 0.7) to the lowest metallic- 
ity ([Fe/H]=— 2.2) isochrone. As stars with intermediate 
metallicities will be interpolated between the template 
populations, this translates into distance uncertainties 
of at most 10%. The other contributions to the distance 
uncertainty are inaccuracies in the isochrones and the 
photometric uncertainties, both of which are expected 
to be much smaller. The overall distance uncertainty is 
therefore of order 10%. 

We deliberately choose to avoid fitting the thin disk, 
because a template population with more typical thin- 
disk populations should have a large range in age and 
metallicity and would be difficult to distinguish from 
a combination of the three templates described above. 
Hence, when presenting and interpreting our results, we 
will ignore the distance ranges where the thin disk is 
expected to be an important contribution to the star 
counts, and only consider regions with Z >1.0 kpc, or 
roughly four thin-disk scale heights. 

All stars with 15< g <21.5, 15< r <21.0 and 0.1< g- 
r <0.8 are used in the fits, with Hess diagram pixel sizes 
of 0.07 in color and 0.2 in magnitude. This resolution 
is sufficient to resolve the features on which the anal- 
ysis depends, but still ensures sufficient signal-to-noise 
per pixel. The color-cut of g-r <0.8 eliminates the is- 
sue of fitting the faint, red thin-disk stars. For each of 
the templates, model CMDs are created for a range of 
distance moduli. The magnitude limits used in the fits 
(15< g <21.5) correspond to distance limit for MSTO 
stars of roughly 1.5< D <25 kpc. Since the MSTO stars 
provide the best constraints for these CMD fits, only this 
distance range will be analysed. However, to avoid edge 
effects, model templates are created for distance mod- 
uli between 5.0 (100 pc) to 20.0 (100 kpc) in steps of 0.2 
mag (~10% distance bins). MATCH then determines the 
best-fitting combination of model CMDs. Subsequently, 
the uncertainties on the obtained results are determined 
by refitting Monte Carlo realizations of e ach CMD draw n 
from its best-fit model, as described in iDolphinl (|2001l) . 
It is again important to note the difference between this 
approach and using stars of a particular spec tral type to 
make 3D maps of the Milky Way (as in e.g. I Juric et all 
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12008ft . The CMD fitting will give a direct, high signal-to- 
noise map of the stellar mass density, even if the stellar 
populations are spatially varying. This approach also uti- 
lizes the information for stars of different spectral types 
(colors, luminosities) simultaneously. 

A comparison of the best-fit model with the CMD data 
for the region at 1=94° and 6=+30° (the same field as 
shown the lower left panel in Figure [2]) is presented in 
Figure 2] Although the model is not perfect (which per- 
haps should not be expected due to the simple model 
population make-up) , the general features of the density 
distribution are well-reproduced. The mean residuals are 
25% of the mean star counts across the portion of the 
CMD used in the fit. 

4. RESULTS 

Rather than the stellar mass or the number of stars, 
the CMD fitting software provides the star formation 
rate (SFR) in M Q yr _1 corresponding to each population 
template for each distance modulus bin. However, these 
SFRs can easily be converted to the spatial stellar mass 
density using 

P *' p = ^3-^3) x SFR, At, (1) 

where SFR p is the star formation rate in M Q yr^ 1 of 
template population p, At is the width of the age bin, 
oj is the solid angle in steradian corresponding to each 
CMD, and D\ and Z?2 are the distances corresponding 
to the limits of the distance modulus bin. The 2001 
MATCH version uses a single power-law initial mass 
functi on (IMF) . In t his analysis we opt for a Salpeter 
IMF (|Salpeter IH955T ) with cutoffs at 0.1 and 120 M , 
which yields a total mass of about a factor 2 high er than 
currently favored IM Fs (e.g. iKroupa et al.lll993f) for old 
populations (see also iMartin et al.l 120081 ) . It should also 
be kept in mind that these masses correspond to the 
total mass of stars formed, which, especially for old pop- 
ulations, does not correspond to the total mass in stars 
today, but rather to the total mass in both stars and 
stellar remnants. 

Figure [3] shows the stellar mass densities obtained from 
the CMD fits for each of the stripes, summed over all 
three metallicities. Figure [6] shows the same, but with 
a color-coding that represents the relative contributions 
of the three different template populations. Generally 
speaking, the stellar density decreases in a continuous 
fashion with increasing distance from the Galactic center 
and plane, although some deviations are noticeable. For 
example, a region of apparent excess density extends to 
the bottom of the ?= 94° panel in Figure El c aused by the 
Sagitt arius stream (|Maiewski et al.1 l2003t lYannv et al.1 
2009b). The contributions of the three different tem- 
plate populations change as expected. With increasing 
distance from the Sun, the 'thick disk', 'inner halo' and 
'outer halo' populations dominate in turn. The Galacto- 
ccntric distance where the metal-poor 'outer halo' starts 
dominating over the more metal-rich inner halo varies 
between 10 and 20 kpc. 

Figure [7] shows the mean density and mass-weighted 
metallicity as a function of distance from the Sun for a 
set of different bins in Galactic latitude. Apart from the 
total mass density, the contributions from the individ- 
ual template populations are also plotted. The increas- 



ing contribution of the thick disk to the density along 
lines of sight with decreasing latitude is clearly visible, 
both in the density graphs and in the steepness of the 
metallicity gradient at distances less than 10 kpc. The 
metallicity graphs indicate that the stellar halo has an av- 
erage metallicity of [Fe/H]~ —1.6 a t distances less tha n 
~15 kpc, confirming the findings of llvezic et al.l (2008), 
but then drops to [Fe/H]< —2 further out. Although 
the transition is smooth, and both 'halo-like' template 
populations contribute to the mass density at every dis- 
tance, such a drop is consistent with th e presence of 
two d istinct populations, as predicted by ICarollo et al.1 
(|2007| ). Recall, however, that the transitions between 
different templates can only be taken as direct evidence 
of a metallicity gradient in the halo, if age differences 
between halo stars play no role in their turn-off color. 
For the thick disk and stellar halo this is a reasonable as- 
sumption, as both are generally believed to be relatively 
old (>10 Gyr), so that the isochrones are much less sen- 
sitive to age differences than to metallicity differences. 

The stellar density values, p* tP , are also presented in 
Table [TJ As the table contains 25 476 rows, only the first 
five rows are printed, with the complete table available 
in electronic form online. 

4.1. Smooth Galactic model fits 

To describe the total stellar densities (i.e. summed 
over all three template populations) quantitatively, we 
fit a model to all data values p*^ represented by the 
maps in Figure [5J Any substructures that are present 
on top of a 'smooth' stellar distribution should become 
apparent after subtraction of the best-fit model. The 
model we use consists of a double exponential thin disk 
and t hick disk (|Bahcall fc Soneirall 19801 : iGilmore fc Reid l 
Il983fl . and an axisymmetric power-law halo (lEggen et al.1 
[19621 iChiba fc Beersll2000t I.Turic et all 121)03 ): 



p(R,Z) 




Here, Rq is the distance from the Sun to the Galac- 
tic center; L thin , H thm , L thick and H thick are the scale 
lengths and heights of the thin and thick disk; qualo and 
n haio are the halo flattening and power-law exponent; 
Pthin,Q is the local thin-disk density and fthick,Q and 
fhaio,Q are the local density fraction of the thick disk 
and halo relativ e to the th i n disk , respectively. In light 
of the results of iBell et al.l (|2008l ) there seems to be no 
need for a more sophisticated, for example triaxial, halo 
model. For a grid of parameter values, models are com- 
pared to the density val ues resulting from the CMD fits. 
Assuming i? Q =7.6 kpc (|Valleej l2008) 7 . the distance and 
the Galactic latitude and longitude of each bin can be 
converted to R and Z. Since the density in each bin, 

7 As Rq is still not very well determined (e.g. Bow et al. 2009), 
we have also used values of 7.1 and 8.1 kpc for Rq and found the 
effects of these changes on the fit results to be negligible compared 
to the uncertainties quoted in Table [3] 
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Fig. 5. — Tomography of the Milky Way. Shown are maps of the stellar mass density in a set of slices of constant I, drawing on the 
single population fits, as function of distance from the Sun, and height above or below the Galactic plane, in kpc. Contours correspond 
to steps of a factor 10 in stellar mass density. Each panel shows a different imaging scan, the Galactic longitude of which is listed at the 
bottom, and darker gray levels correspond to higher densities. The dashed semicircles show the maximal distance to which fit results are 
(mainly) based on main-sequence turn-off star colors and densities. 
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Fig. 6. — Density-metallicity tomography: shown is the distribution of stellar populations from the triple population fits with the 
same layout as in Fig. [5] The thick-disk-like population with [Fe/H]~ —0.7 is color-coded as red, the inner halo-like population with 
[Fc/H]~ —1.3 as blue, and the outer halo-like population with [Fc/H]~ —2.2 as green. The white, dashed semicircles show the distance up 
to where fit results are (mainly) based on main-sequence turn-off star colors and densities. 
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TABLE 1 

Stellar density maps 



I b m-M pi o-pi p 2 u p2 p3 <Tp 3 



(°) 


(°) 


(mag) 


(M Q kpc- 3 ) 


(M Q kpc" 3 ) 


(M Q kpc" 3 ) 


(M Q kpc- 3 ) 


(M Q kpc" 3 ) 


(M Q kpc" 3 ) 


110.0 


-34.5 


11.1 


0.0 


4.76E5 


0.0 


5.00E3 


1.28E7 


5.12E6 


110.0 


-34.5 


11.3 


0.0 


1.16E5 


0.0 


6.86E5 


3.14E6 


3.01E6 


110.0 


-34.5 


11.5 


0.0 


1.51E5 


0.0 


2.75E5 


5.74E6 


1.87E6 


110.0 


-34.5 


11.7 


1.11E5 


1.34E5 


0.0 


7.00E4 


0.0 


1.28E6 


110.0 


-34.5 


11.9 


2.36E5 


2.00E5 


0.0 


6.49E4 


5.1E5 


6.46E5 



Note. — Densities pi, p2 and P3 correspond to the mass densities for the template populatio ns with [Fe/ H]=— 2.2, 
—1.3 and —0.7, respectively. The densities in this table correspond to the Salpeter IMF ((Halpctcr 1955]) used by 
the MATCH software and are in solar masses per cubic kilopars ec, unlike the densities quoted in the text of the 
paper, which are corrected for a Kroupa IMF (Kroupa ct al. 1993) and in solar masses per cubic parsec. This table 
is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance 
regarding its form and content 
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Fig. 7. — Mean stellar mass densities and metallicities as function of distance. The left-hand panels show the total stellar mass density 
with the black error bars for four different bins in Galactic latitude, indicated in the top right of each panel, and averaged over all longitudes. 
Colored lines show the mass density in the individual template populations with red for the 'thick-disk-like', blue for the 'inner-halo-likc' 
and green for the 'outer-halo-like' population. In the middle panels the colored lines indicate the fractional contribution of the individual 
template populations to the total mass density on a scale of to 1, following the same color scheme as in the left-hand panels. The 
right-hand panels show the mass-weighted mean metallicity. 
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p*^i(R, Z) has an uncertainty from the Monte Carlo tests, 
the x 2 goodness-of-fit of each model can be determined 
straightforwardly as 



X 



E 



pi 



odel 



Pi,obs 



(3) 

Oi observed ' 
l 

As mentioned before, we only use the bins for which 1.5< 
D <25 kpc and Z >1.0 kpc. 

Some localized stellar overdensities are known to ex- 
ist in the Milky Way which may be strong enough to 
affect significantly the attempt to fit a 'smooth' model. 
The strongest one in the halo is the Sagittarius stream, 
presumed to be tidal debris from the Sagittarius dwarf 
galaxy. It arcs through the North Galactic cap through 
the str ipes at /=229°, 270° 300°, and 330° at high lati- 
tudes (|Belokurov et al.l [20061 ). and inter sects the stripes 
at ^==94° and 178° at nega tive latitudes (jMajewski et al.1 
l200l lYannvlt~aIll2009bl) . This can already be seen in 
the density maps in Figure EJ for example at Z=94°, 
Z < —15 kpc and D <10 kpc. At lower latitudes the 
largest known overdensity is the Mo noceros overdensity, 
discovered by iNewberg et al.l (|2002D and of still contro- 
versial origin. It may plausibly be part of the Low Lati- 
tude or anti-center stream, a ring- like structure though t 
to encircle the Milky Way completely (jlbata et al.| [2003). 
This structure can also be seen in Figure [5l near the 
position of the original discovery (in particular for the 
;=178° and /=187° stripes, where a 'bump' is apparent 
in the stellar density at 5< D <10 kpc and 0< Z <5 
kpc). 

To test for the influence of these structures on our fit 
results, we fit the density maps twice, once with all data, 
and once without the regions dominated by these struc- 
tures. The regions that are excised to remove the impact 
of the Sagittarius stream and Monoceros are listed in 
Table [2] amounting to of all data points. 

The role of the thin disk is limited by only looking at 
the area at \Z\ >1.0 kpc, and we fix its pa rameters at 
Lthin= 2-6 kpc and H t hin=0-25 kpc, following I Juric et al.l 
(12008T) and in lin e with earlier star count results (e.g 
ISiegel et alJl2002t ). In other words, we treat the thin 
disk as a known and fixed 'perturbation' to our CMD 
analysis. We are now left with seven free parameters in 
fitting the total stellar mass density (summed over the 
three metallicity components): the thick-disk parameters 
f thick, Lthick and H thickl the halo parameters f ha i a , 

Qhalo 

and nhaioi and the local thin-disk density Pthin,s- The 
latter serves as normalization for the thick disk and halo 
density, hence in practice reducing the number of free 
parameters to six. 

4.2. Structural parameters of the thick disk and stellar 

halo 

In Table [3] we list the parameter values correspond- 
ing to the best fit s, with the mass de nsities scaled to a 
Kroupa-like IMF (jKroupa et al.l 1 1993ft . The best fit to 
all data has a reduced \ 2 °f 4.2, while the best fit to 
the data without Sagittarius and Monoceros has a re- 
duced x 2 of 3.9. Inclusion of these structures thus has 
a noticeable impact on the goodness-of-fit, but in both 
cases the smooth mo dels provid e a po or de scription of 
the da ta, in line with iBell et all (|2008| ) and I Juric et "all 
(2008). Since this inhibits deriving proper confidence in- 
tervals based on the x 2 statistic, the uncertainties on the 
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Fig. 8. — Results of fitting smooth Galactic models to the stellar 
density maps. For the case where all data are used, the best fit 
values are indicated with squares and iso-x 2 contours are plotted 
with solid lines. Triangles and dashed contours are for the case 
when the Sagittarius stream and the Monoceros overdensity are 
excised. In both cases, contours are scaled to correspond to the 1, 
2 and 3<r uncertainties determined through bootstrap tests. 
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Fig. 9. — Thick disk (top) and halo (bottom) densities in the 
best-fit smooth model as function of distance from the Galactic 
plane. The densities were calculated for two galactocentric radii, 
7.6 kpc (the assumed galactocentric radius of the Sun) and 12.5 
kpc, and the sequences of points are labeled accordingly. Dots 
show the density at each location for the best-fit model, while the 
grey areas show the 95% confidence region determined from the 
bootstrap analysis. 



fit parameters listed in Table [3] were obtained by resam- 
pling the data points, p* t i(R,Z), as in bootstrap tests. 
For each test, resampled data sets were created by ran- 
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domly drawing data points (with replacement) from our 
set of stellar densities, p*,i{R, Z), until the same number 
of data points is reached. The fit parameters were then 
re-determined for the resampled data sets by re-fitting 
the model given by Equation [2] This procedure was re- 
peated 200 times for the full data set and for the data 
set without Sagittarius and Monoceros. In Figure [5] we 
present the contour plots of the la, 2cr, and 3a (68%, 
95% and 99.7%) confidence levels around the fit param- 
eters, which are also based on the bootstrap analysis. 

The mass densities of the thick disk and halo are very 
well constrained, as demonstrated in Figure [9j where 
the best model densities and their 95% confidence in- 
tervals are plotted at different locations. The thick-disk 
density is best constrained at (Rqc ,\Z\)=(7.6,2.5) kpc 
as pthick = 10~ 3 - 70±0 - 02 M pc~ 3 , and the halo density 
at (i? GC ,|Z|)=(10.0,12.5) kpc as p halo = lO" 519 * 001 
M pc -3 (i(r 5 - 21±0 - 01 M pc~ 3 when Sagittarius and 
Monoceros are omitted). These formal errors are pre- 
sumably exceeded by uncertainties in the IMF and as- 
sociated with substructure. At the solar location the 
uncertainties on the density are significantly larger than 
at the locations stated above, as the local values are not 
directly constrained by the data. Rather, they are ex- 
trapolations from data at higher \Z\ and Rgc, and are 
therefore model dependent. However, since most previ- 
ous work has focussed on local samples of stars, we will 
use the local extrapolations of our results in order to com- 
pare with earlier studies. The value we find for the local 
mass density of the thick disk is pthick.Q = 10~ 2 - 3±OA 
Af Q pc- 3 , and for the stellar halo p ha io,Q = l()- 4 - 20±0 - 05 
M Q pc -3 (pmck,e = 10- 2 - 3±0 - 1 M pc -3 and p h alo,G = 
-LQ-4.i8io.05 jVf Q p C -3 wnen Sagittarius and Monoceros 
are omitted). 

Recent dynamical modeling (e.g. lFlvnn et al.ll2006D im- 
plies a local thick-disk density of lO _2 46 M pc -3 and a 
local halo density of lO _4 M pc -3 , similar to our find- 
ings. However, our local mass densities include stellar 
remnants, and a proper comparison should take this into 
account. The total local mass density, dominated by the 
thin disk (in stars and white dwarfs), is ~ O.O38Af pc -3 
(|Flvnn et al.l [2006) . Star-count studies typically yield 
thick-disk and halo normalizations with respect to the 
thin disk of 1 to 10% and 0.1 to 0.2%, respectively 
(|Siegel et al.l 120021) . This therefore implies a local den- 
sity, including remnants, of 10~ 3 - 40 to 10 -2 42 M pc -3 
for the thick disk and 10" 442 to lO" 412 Afopc" 3 for the 
stellar halo, values that are consistent with our density 
estim ates. If we normaliz e our findings by the local value 
from iFlvnn et al.l (120061) . we find local thick-disk and 
halo mass-d ensity fractions of 15±4% and 0.17±0.03%, 
respectively. Uuric et all (|2008D . in their analysis of SDSS 
data, find 12±1% for the thick disk, consistent with our 
results, and 0.51±0.13% for the halo, inconsistent at the 
2.6<7 level. These normalizations also imply that the stel- 
lar density of the thin disk should never amount to more 
than ~10% at \Z\ >1.5 kpc in our maps, confirming that 
it is merely a perturbation to our analysis at these loca- 
tions. 

Thick-disk structural parameters determined previ- 
ously throu gh star-count stu dies have shown a large 
spread (e.g. ISiegel et al.|[2002D . Scale heights vary from 
^0.6 to ~2.0 kpc, with the more recent measurements 



TABLE 2 

Regions Excised for Second Model Fit 



Structure 


I (°) 


b(°) 




m — M (mag) 


Sgr stream 


94 


< -65 




>16.0 


Sgr stream 


178 


-65 to - 


25 


>16.5 


Sgr stream 


229 


+55 to 4 


70 


>15.5 


Sgr stream 


270 


+55 to 4 


70 


>15.5 


Sgr stream 


300 


+58 to 4 


72 


>15.5 


Sgr stream 


330 


460 to 4 


80 


>15.5 


Monoceros 


178 


+15 to 4 


35 


14.0-16.0 


Monoceros 


187 


415 to 4 


35 


14.0-16.0 



converging to the lower end of this range; and scale 
lengths vary from 2.8 to 4.3 kpc. Uuric et all (|2008f) find 
very different values depending on the tracer stars or 
photometric parallax relations that are used, but settle 
on 0.9±0.2 kpc and 3.6±0.7 kpc for the thick-disk scale 
height and length. Through our present analysis we have 
been able to estimate the scale height and length of the 
thick disk to 10% accuracy: H t hick — 0.75 ± 0.07 kpc 
and Lthick = 4.1 ± 0.4 kpc. This is consistent with but 
more precise than previous determinati o ns, an d in good 
agreement with the work of Uuric et al.l (|2008f) . 

Halo parameters, when described by a sin gle power- 
law, are probably varying with radius (e. g. iBell et all 
120081: IWatkins et al.|[2009t ISesar et al.ll2010D . and in this 
study only the stellar halo within a radius of 25 kpc from 
the Sun has been considered. Two recent studies based 
on SDSS data have measure d the flattening a nd power- 
law index of the stellar halo. Uuric et "ail (|2008l ). reaching 
to ~20 kpc, find a halo flattening parameter of %=0.64 
and a power-law index of rih=— 2.77 w ith quoted uncer- 
tainties of <0.1 and <0.2, respectively. IBell et al.l (120081 ) 
find that a qu of 0.6 to 0.7 gives the least excess rms scat- 
ter around their model fits, but that a qt of 0.7 to 0.8 
gives the best x 2 values. In either case, the constra i nts on 
qu are weak. The halo model used by I Bell et al.l (2008) 
is a broken power-law, with independent indices on ei- 
ther side of the break radius. These parameters are also 
weakly constrained, with the preferred break radius ly- 
ing between 20 and 30 kpc, the inner power-law index 
between 2 and 3, an d the outer power-law index between 
3 and 4. Similarly, IWatkins et ail (j2009[ ) use RR Lyrae 
variable stars in SDSS Stripe 82 to infer rih=— 2.4 within 
23 kpc and rih—— 4.5 at larger distances. By compar- 
ison, our analysis yields very precise measurements of 
n h = -2.75 ± 0.07 and q h = 0.88 ± 0.03 for the stel- 
lar halo within 25 kpc from the Sun. While our value 
for the power-law index is in excellent agreement with 
Juric et all (120 081 and the inner power-law index from 
Bell et all (120081 ). we find the stellar halo to be less flat- 
tened. Note that these previous SDSS determinations 
were restricted to b >30°; the inclusion of the SEGUE 
scans across the Galactic plane boosts the power to de- 
termine flattening. 

4.3. Substructure, or deviations from a smooth model 

To focus on substructure in the stellar halo and thick 
disk we now subtract the best-fit smooth model to the 
complete data set from the density maps in Figure[5] The 
residuals after subtracting this smooth model are shown 
in Figure fTOl revealing a wealth of substructure. In par- 
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TABLE 3 
Smooth Galactic model fits 



Parameter 


Full Data 


Exclude Sgr.Mon 


Pthick,Q (M e pc- S ) 


1Q-2.3±0.1 


ig-2.3±0.l 


Lthick (kpc) 


4.1±0.4 


4.0±0.4 


Hthick (kpc) 


0.75±0.07 


0.75±0.08 


Phalo,Q (M pc- 3 ) 


2Q-4.20±0.05 


1 g-4.18±0.05 


Qhalo 


0.88±0.04 


0.85±0.03 


n halo 


-2.75±0.07 


-2.80±0.07 



ticular, many overdense regions jump out, and although 
there are also underdense regions, the density contrast of 
the latter with respect to the smooth model appears to be 
much weaker. In principle the CMD fitting method can 
provide not only the stellar densities, but also yield infor- 
mation on the stellar populations of the overdensities, for 
example their metallicities. However, the current analy- 
sis uses three simple template populations, which differ 
in metallicity but all have the same age (~14 Gyr). For 
the thick disk and stellar halo it is justifiable to interpret 
the population differences detected in our fits as metal- 
licity differences, as both components consist of old stars. 
Overdensities, however, might be streams populated by 
extra-galactic sources or material swept up from even the 
thin disk. There is no reason to believe that their stars 
should have the same ages as the smooth components 
on which they are superposed. The reconstruction of 
the detailed stellar population properties of specific sub- 
structures thus requires a more detailed analysis of each 
overdensity. Therefore, we limit ourselves here to listing 
the overdensities we detect and a short description of the 
main structures responsible for the overdensities seen in 
Figure Q33J 

As the strength of an overdensity detection depends 
on the smooth model that is subtracted from the density 
maps and each model parameter has a specific uncer- 
tainty, we need to define what we consider a significant 
detection. We do this by constructing a set of smooth 
Galactic models on a grid of parameters, and selecting 
all models which are within ler of the best-fit model in 
the 6D parameter space given by the model parameters 
in Table El Density maps are evaluated for all pixels in 
these 282 models. The pixels with densities at least 3ct 
higher than the model density for more than half (141) 
of the models are considered to be significant. Table 
[4] lists all significantly detected overdensities, most of 
which consist of a large number of such pixels. The peak 
significance listed in the table corresponds to the most 
significant pixel in each overdensity with respect to the 
best-fit smooth model. The total significance of the over- 
density is much larger in most cases, as they consist of 
many pixels. Because the metallicities and ages of the 
stars in the overdensities do not necessarily fall neatly 
within one of the template populations used in the CMD 
fits, the distances listed should be taken as estimates, 
not as precise measurements. Note that one detection is 
caused by two bright globular clusters, NGC 5024 (M53) 
and NGC 5053, which happen to lie exactly in the scan 
at Z=330°. 

Sagittarius stream. Probably the most striking stel- 
lar 'substructure' in the Galactic halo, the Sagittarius 
stream (labeled with an 'S' in Fig. H0| is visible in several 



stripes, both in the northern and the southern Galactic 
hemisphere. The northern arm is clearly visible in the 
stripes at Z=229°, 270°, and 300°, just outside the dashed 
semicircles marking a dis tance of 25 kpc. At 1= 330° the 
stream is more distant (jBelokurov et alJ l2006| ) and off 
the plot . The two arms in wh ich this part of the stream 
is split (|Belokurov et alJfeOOGt ) are detected very clearly, 
particularly in the Z=270°stripe. Typical peak densities 
in these stripes are 5 to 8xlQ~ 6 M pc -3 . At Z=203° the 
northern Sgr stream should en ter the stripe around 
6=30° and plunge into the disk (|Belokurov et alJ [20061 : 
lYannv et~a l. 2009b) . As it is roughly equidistant with the 
Low-Latitude Stream, it is difficult to distinguish it in the 
overdensity at low latitudes in this stripe. In the south- 
ern Galactic cap, the Sgr stream yields strong detections 
in two stripes, namely the 1=94° and 1 7 8° strip es, re- 
producing the detection by lYannv et alJ (|2009bt ). Here 
the stream is also seen at distances of 20 kpc and higher. 
The density of the stream in these southern detections 
peaks at ~lxlO~ 5 M Q pc _3 , slightly higher than in the 
North Galactic cap. 

Virgo. The Virgo overdensity (|Juric et al.ll2008D is an- 
other strong overdensity in the Northern Galactic cap, 
centered at (?,&)~(280°,70°) and extending between dis- 
tance s of -10 to -20 kpc (jBell et alj[2008t Uuric et all 
2008). It is very prominent in the residual maps of the 
stripes at ^=270° and 300° in the same direction as the 
Sgr stream, but at distances less than 20 kpc. The 
overdensities at similar distances in the /=330° stripe 
might be related to Virgo as well. Typical densities 
in the Virgo overdensity we find in the ^=270° and 
300° stripes are ~2xlO _5 M pc -3 , peaking at around 
~4xlO~ 5 Ai0 P c~ 3 . 

Hercules-Aquila. A large overdensity extending both 
above and below the disk , and spanning —80° in lon- 
gitude, was detected by iBelokurov et al.l (|2007l ) and 
christened the Hercules-Aquila cloud. This struc- 
ture is centered at I ^40° but is visible out to 
I <90° (jBelokurov et al.l 120071) . The overdensities de- 
tected in the stripe at Z=70° and distances larger than 
—22 kpc might therefore be associated with this recently 
discovered stellar component. 

Monoceros and Low Latitude Stream. The Mono- 
ceros structure was first discovered in SDSS data by 
iNewberg et alJ (|2002j ) as a strong stellar overdensity at 
(I, b) —(200°, 20°). Later such overdensities were detected 
over a large range of Galactic longitudes both above and 
below the Galactic plane, always at low latitudes and 
roughly equidistant, thu s forming a ring- l ike structure 
around the Milky Way (llbata et all 120031 : lYannv et all 
120031 : IConn et all 120051 120071 ). Ta the stripes at the 
longitudes close to the original discovery, Z=178°, 187°, 
203°and 229°, the overdensity is detected at high signif- 
icance at heliocentric distances of roughly 10 kpc and 
—3 kpc above the Galactic plane. The ring-like nature 
of the Low Latitude Stream appears to be confirmed 
by our new map, as overdensities at the same helio- 
centric distances and latitudes are present in the stripes 
at ;=94°, 130° and 150°, but not at Z=110°. The den- 
sity of the structure is highest in the stripes at Z=178°, 
187°, and 203°, where the maximum density reached is 
— 1.5xlO~ 4 Af pc~ 3 . The density appears to fall off by 
a factor of —2 in the stripes at both lower and higher 
longitudes. The overdensity detections at negative lati- 
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Fig. 10. — Residual stellar mass density after subtraction of a smo c model. Panel layout is the same as in Figure [5] Starting 

from black, the grayscale levels correspond to areas with residual densities > 3, 3—2, 2—1, 1 to — 1, and < —1 times the model density. 
White areas contain no information due to an absence of data or too high reddening. The dashed semicircles show the distance up to where 
fit results are (mainly) based on main-sequence turn-off star colors and densities. Four known overdensities are labeled: the Sagittarius 
stream (S), the Virgo overdensity (V), the Hercules- Aquila cloud (H), and the Monoceros, or Low Latitude stream (M). The panels for 
2=270°, 300°and 330° have insets showing the inner 5x5 kpc region. 



tudes (jlbata et al. 2003: lConn et~aT1 l2005. 2007) are con- 
firmed here, with an overdensity detected consistently 
-4 kpc below the plane from Z=130° to ^203°. Be- 
low the plane the overdensities have lower densities than 
above the plane, increasing from ^O.2xlO -4 M0pc~ 3 at 
Z=130° to ~1.3xlO~ 4 M^pc- 3 at 1 =130°. 

Detection at (R,Z)=(6.5,1.5). Uuric et all (|2008l ) 
detected an overdensity at Galactocentric coordinates 
(R,Z)=(6.5,1.5) kpc. In the stripes at Z=270°, 300° and 
330° we find a very nearby overdensity at b ~65° and at 
a distance range from 4 kpc to the smallest distance we 
probe, 1.5 kpc. This over density therefore c orresponds to 
the same as detected by Uuric et al.l (|2008|) . The nature 
of this overdensity is so far unknown. 

Unknown substructure. In a few locations less promi- 
nent substructures are visible, that cannot directly be 
identified with known overdensities. Most of these are 
located at low latitudes and might be related to the Low 
Latitude Stream. In a future publication we will study 
these overdensities and their possible connection with the 
Low Latitude Stream in detail. 

5. SUMMARY AND CONCLUSIONS 

We have applied CMD-fitting techniques to SEGUE 
photometric data to study the thick disk and stellar halo 
of the Galaxy at both high and low latitudes. Three 
template stellar populations, all with an age of ~14 Gyr 
and metallicities of [Fe/H]=— 0.7, —1.3 and —2.2, are fit 
to the data, yielding a three-dimensional census of stel- 
lar mass. Since the thick disk and stellar halo consist 
of predominantly very old stars, differences in turn-off 



color are indicative of differences in metallicity. A change 
of turn-off color is seen in the halo at distances of 15 
to 20 kpc, which is therefore interpreted as a change in 
metallicity. At distances smaller than ~10 kpc our re- 
sults indicate a mean metalli c ity of [Fe/H]~ —1.6, similar 
to the work by llvezic et al.l (|2008t ) , who find a value of 
[Fe/H]~ —1.5 out to ~9 kpc. At larger distances, D>15 
kpc, our results indicate that the mean metallicity of the 
stellar halo is [Fe/H]~ —2.2, however. This transition is 
consistent wit h the inference from local samples of stars 
(D<4 kpc~) bvlCarollo et ail J2007T) and preliminary re- 
sults from iBeers et al.l ( 2009T ) , but this is the first time 
this change in population properties is measured in situ. 

Structural parameters of the Galaxy can be derived 
from the resulting stellar-density maps through the com- 
parison with models. Our fits of models with an expo- 
nential thin and thick disk and a power-law halo (Eq. [2]) 
yield constraints on the thick disk local density, pthick,Q, 
scale height, H t hick, and scale length, L t hick, and the halo 
local density, Phaio,Q, flattening, qhaio and power - law in - 
dex, nhaiq . As point ed out by iBelokurov et al.l (|20Q6Ik 
IBell et al.1 d200l . and Uuric et all (l2008h . the SDSS data 
of the stellar halo contains clear evidence of the pres- 
ence of_suhs^ucture with respect to smooth halo mod- 
els. IBell et all (|2008l ) also demonstrated that the larger 
known substructures in the halo significantly influence 
their model fits. To test the sensitivity of our model fits 
to known substructure, we perform our fits to the full 
data set, as well as to a data set with the Sagittarius 
stream and Monoceros overdensity removed. The results 
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TABLE 4 
Detected overdensities 



I (°) 




D (kpc) 




ID 


i 


( — 23,-21) 


(14,16) 


5. / 




70 


(+36, +38) 


(15) 


8.5 




70 


(—23,-22) 


(22,25] 


1.8 


Hcrculcs-Aquila? 


70 


( — 13,-10] 


(23,25] 


3.1 


Hercules- Aquila? 


70 


[+10, +12) 


(23,25] 


3.5 


Hercules- Aquila? 


9-1 


(—82,-70) 


(21,25] 


6.7 


Sgr stream 


94 


(+18, +22) 


(12,14) 


4.1 




91 


(+32, +35) 


(14,16) 


8.4 




130 


(—18,-16) 


(13) 


1.2 




130 


(+19, +22) 


(8,13) 


4.7 




130 


(+26, +27) 


(9) 


7.1 




130 


(+35, +37) 


(9,13) 


5.0 




150 


(—31,-21) 


(11,14) 


3.9 




150 


f i i T i on \ 

(+17, +29) 


(8,13) 


1 10 


Monoccros 


178 


(—62,-39) 


(22,25] 


6.2 


Sgr stream 


178 


(—59.5) 


(14) 


3.1 




178 


(—48,-46) 


(11.5) 


3.5 




178 


(—42,-40) 


(12.5) 


3.2 




178 


(—30,-26] 


/I 1 1 K\ 

(11,15) 


8.3 




178 


(+16, +32) 


(5,14) 


9.7 


Monoceros 


187 


(—30,-17) 


(8,13) 


5.7 




187 


(+14, +35) 


(6,15) 


13.1 


Monoceros 


203 


(—28,-15] 


(8.5,13) 


8.5 


Monoceros 


203 


(+14,+3l) 


(7,15) 


22 


Monoceros 


9n^ 
zuo 






o . o 




229 


(+15,+40) 


(7,18) 


7.1 


Monoccros 


229 


(+35, +38) 


(21,24) 


6.5 




229 


(+57,+64) 


(23,25] 


4.2 


Sgr stream 


270 


(+58,+69) 


(23,25] 


4.1 


Sgr stream 


270 


(+55,+74) 


(2,23) 


8.8 


Virgo/(6.5,1.5) 


300 


(+65, +70) 


[1-5,4) 


5.4 


(6.5,1.5) 


300 


(+58,+72) 


(10,19) 


16.1 


Virgo 


300 


(+60,+65) 


(21,25] 


4.1 


Virgo/Sgr stream 


330 


(+57,+72) 


[1.5,4) 


5.5 


(6.5,1.5) 


330 


(+55,+73) 


(12,24) 


7.2 


Virgo 


330 


(+78.+81) 


(16.5,22) 


12.2 


NGC5024/NGC5053 



Note. — For Galactic latitude b and distance D ranges are 
given within brackets, with square brackets indicating an edge 
of the data rather than the edge of a detection. 



are tabulated in Tabled and are not strongly influenced 
by the presence of these structures. In either case, the 
reduced x 2 of the best fit, 4.2 and 3.9, shows that the 
smooth model is a poor fit to the data. 

The values we find for the local density of the thick 
disk, pthick,Q = 10~ 23±01 M Q pc -3 and the local halo 
density, phaio,Q — 10~ 4 - 20±005 M pc~ 3 , agree well with 
star-count studies. In addition, our results for the thick- 
disk structural parameters, H t hick = 0.75 ± 0.07 kpc 
and Lthick = 4.1 ± 0.4, are well with in the range of pa- 
rameters found in previous studies (Sicgcl e t al.l I2002T) 
and consistent with the results of lJuric et all (|2008f ). 
For the flattening of the stellar halo, where we find 
qhaio — 0.88 ± 0.04, our analy sis see ms to diver g e from 
the work by iBell et al.l (|2008l ) and lJuric et al.l (|2008h . 
who use SDSS data to study the halo out to similar dis- 
tances. On the other hand, the power-law index we re- 
cover, rihaio = —2.75 ± 0.07, is in excellent agreement. 
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Subtracting the best model from the stellar-density 
maps unveils abundant substructure (Fig. [TUJ) - Most of 
the overdensities seen can be attributed to know struc- 
tures in the halo and outer disk of the Galaxy: the 
Sagittarius stream, Virgo overdensity, and the Mono- 
ceros overdensity or Low Latitude Stream. Whereas the 
origin of the Sagittarius stream is known to be the result 
of the disruption of the Sagittarius dwarf galaxy, the na- 
ture of the other two entities is currently under debate, in 
particular the interpretation of the overdensities at low 
Galactic latitudes remains contentious. A detailed anal- 
ysis of all substructures is outside the scope of this paper, 
but an in-depth study of the low latitude substructure is 
planned to be presented in a future publication. 
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